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ABSTRACT 

We present H I observations of the edge-on galaxy NGC 891. These are among the deepest 
ever performed on an external galaxy. They reveal a huge gaseous halo, much more extended 
than seen previously and containing almost 30% of the H I. This H I halo shows structures on 
various scales. On one side, there is a filament extending (in projection) up to 22 kpc vertically 
from the disk. Small {Mm > \(f Mq) halo clouds, some with forbidden (apparently counter- 
rotating) velocities, are also detected. The overall kinematics of the halo gas is characterized 
by differential rotation lagging with respect to that of the disk. The lag, more pronounced at 
small radii, increases with height from the plane. There is evidence that a significant fraction of 
the halo is due to a galactic fountain. Accretion from intergalactic space may also play a role 
in building up the halo and providing low angular momentum material needed to account for 
the observed rotation lag. The long H I filament and the counter-rotating clouds may be direct 
evidence of such accretion. 

Subject headings: galaxies: halos — galaxies: kinematics and dynamics — galaxies: structure 
— galaxies: individual (NGC 891) 



1. Introduction 



Recent deep observations of the neutral hydrogen of several nearby spiral galaxies indicate that up to 
about 15% of the neutral hydrogen of a spiral galaxy is located in the halo. The best examples are: NGC 891 
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dSwaters. Sancisi. & van der Hulslll997l): NGC2403 Jschaap. Sancisi & Swatersl200ol : Fratemali et alboOlb : 



NGC4559 dBarbieri et al 



2005b: M31 dwestmeier. Braun. & Th ilkei^lOOs'): UGC7321 (' Matthews & Wood 
2003h: NGC 253 (iBoomsma et all boOSl J): NGC 5775 Lee et all 1200 ID and NGC 6946~ (lBoomsma et al 



2005 al : iBoomsma 1 120070 . Overall, the kinematics of the halo gas is quite regular. The main motion is 
differ ential rotation parallel to the plane. However, the rotation velo cities in the halo are lower than in the 



disk (|S waters. Sancisi. & van der Hulstill997|^ 



Fratemali et al 



2002h . In some cases a small overall radial 



inflow is found superposed on the rotation (IFratemali et al.ll2001). On small scales stron g vertical motions 
from and towards t he disk are also observed (Boomsma et al. 2005a; Boo msma 1 120071 ). Halo gas is also 
detected in Ha (e.g. ^Hoop es. Walterbos & Randil999i:iRossa & Dettmarii2003l). with kinematics similar to 



that of the neutral gas (iRandlbood: 



TuUmann et al 



2006al. 120071 : iKamphuis et al.ll2007h . and in X-rays (e.g. lStrickland et al.ll2004^ . 



200i£ 
(e.g.lsM 



Fratemali. Oosterloo & Sancisill2004l : lHeald et al. 



The origin of the gaseous halos is still a matter of debate. The galactic fountain mechanism (|Shapiro & Field 



19761 ) has received most attention to date. In this scheme, gas is pushed into the halo by stellar winds and 
SN explosions, mostly in t he hot ionized p hase. This gas travels through the halo, eventually cools to neutral 
and falls back to the disk (Bregman 198ol) . There is strong observational evidence supporting the fountain 
mechanism, such as the close correlation between the distribution of Ha and the high- velocity gas found in 
NGC 6946 ( Boomsma II2OO7 ), as well as that the star formation rate appears to correlate with t he luminosity 



of the X-ray halo ( TuUmann et al. 2006bf) and with the radio continuum emission of the halo ( Dahlem et al. 



20061) . However, the sample of galaxies studied so far is still too small to be able to make any statement on 
whether flows related to star formation are the dominant mechanism or not in particular for the neutral halos. 
Furthermore, there are indications th at an ionized halo only builds up if the star formation rate is above a 
critical value jTulhnann et al. 2006b) which may pose a problem for the H I halo d etected in the super- thin 
LSB UGC 7321, a galaxy with a low star formation rate ([Matthews & Woodll2003h . The main problem for 
the fountain mechanism comes, however, from the study of the kinematic s. Simple (i.e. so-called ballistic) 
models have been unable to re produce the kinemat i cs of t he ionized gas dCoUins. Benjamin & Randll2002l : 
Heald et al.ll2006al) . Recently, IFratemali & Binneyl (120061) have shown that also the kinematics of the neu- 
tral halo gas cannot be explained by "pure" ballistic galactic fountains. This suggests that other effects, 
such as the interaction with a pre-existing hot halo or the accretion from intergalactic space, must play an 
important role. There have also been attempts to model th e extra-planar gas as a stationary medium in 
hvdr ostatic equilibrium (iBeniami ni booi iBamabe et al. l2005l) . or as a cooling flow accretion in a CDM con- 
text ( Kaufmann et al hood) , but none of these models is able to reproduce the observations completely. It is 
possible that the halo gas is the result of complex phenomena involving both intemal and extemal processes. 

Understanding the origin and nature of the gaseous halos surrounding spiral galaxies is important for 
several reasons. First, the halo is the region where material can be exchanged between different parts of 
the galaxy and this circulation of gas is fundamental for the galactic hfe-cycle. Secondly, galactic halos are 
the interface between the galaxy, which is visible and well studied, and the Intergalactic Medium (IGM), 
the content and properties of which remain largely unk nown. CDM cosmological models predict tha t most 
of the baryonic material is currently in the IGM (e.g. IWhite & Frenklll99ll : ISommer-LarsenI l2006l ). The 
discovery of the halo gas may provide a new and efficient way to probe the IGM by studying the exchange 
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of material between galaxies and their environment. 



We believe that the halo gas observed in external galaxies is the analogue of the I ntermediate- and High- 
Velocity Clouds (IVCs and HVCs) of the Milky Way dWakker & van Woerdenll 19971) . The cloud complexes 
with anomalous velocities found in galaxies like NGC 2403 ( Fratemah et al.ll2002 ) have similar masses and 
velocity deviations with respec t to the disk as the largest galactic HVCs for which the distances are known 
(e.g. complex A. ,Wakkeiil2001r) . The total gas mass of the HVCs, if located in the halo at distances up to a 
few tens of kpc, would b e of the order of 10^ Mp, similar to tha t of the H I f ound in the halos of e xternal 
galaxies (e.g. NGC 891, Iswaters. Sancisi. & van der Hulstlll997l : NGC2403, Ipratemah et aDl2002l) . It is 
therefore interesting to note th at some of the HV Cs, in particular complex C, have been found to have a low 
metallicity (Z ~ 0.1 - 0.3 Zc^\ iTripp et aLll2003h . This may indicate that some of the HVCs are accretion 
from the surrounding IGM of "unprocessed" materi al onto the disk of our Galaxy. Such accretion may be 
necessary to explain the evolution of disk galaxies ( Naab & O strikei^bood ). Accr etion onto disk galaxies 
may also occur as the merging of small, gas-rich satellites (|yan der Hulst & Sancis 3[l988, 2005). Recently, a 
search has been carried out for the two largest members of the Local Group other than the Milky Way: M3 1 
and M33 (|Westmeier. Braun. & Thilkeii l2005h and a populat i on of high-velocity clouds at large distances 
from these galaxies (about 50 kpc from M31, iThilker et al.l (120041) ) has been found. These clouds have 
typical masses of a few times 10^ Mq. Most of them are thought to be remn ants of the accretion of small 
companion galaxies, although a small fraction might be primordial gas clouds dWestmeier. Braun. & ThiUcer 
2OO5I) . 



In this paper, we present very deep H I observations of the nearby edge-on spiral galaxy NGC 891 
obtained with the upgraded Westerbork Synthesis Radio Telescope (WSRT). NGC 891 is an Sb/SBb galaxy, 
one of the best studied nearby edge-on gala xies. The disk of NGC 891 shows intensive star formation at a 
rate of ~ 3.8 Moyr~' (IPopescu et al.ll2004b . The halo region has b een studied at various wavelengt hs and 
shows a variety of components from radi o continuum emission (e.g. lAUen. Sancisi & Baldwinlll978r) to hot 
diffuse gas (e.g. iBregman & Pildislll994l). NGC 891 is cons idered to be very similar to the Milky Way with 
regard to mass and stellar components (Ivan der Kruitil 19841) although it has a higher star formation rate and 
significantly stronger radio continuum emission. A summary of its physical parameters is given in Tabled] 

In the past, NGC 891 has been studied in H I several times with ev er increasing sensitivity (e.g. 



Sancisi & Allen 



19791 : lRupenlll99ll : ISwaters, Sancisi. & van der Hulstlll997h . The first indication of extra 



planar material had been reported already bv lSancisi & AUenl (Il979h. but thei r favored explanation was that 
of a flaring outer disk. Almost two decades later lBecquaert & CombesI (Il997h proposed a different interpre- 
tation in terms of a warp of the H I disk along the line of sight. Final ly, more sensitive WSRT observations of 
NGC 891 revealed a much more extended extra-planar component ({s waters. Sancisi. & van der Hulst 1997) 
and a careful modeling of the full data cube showed that the most likely explanation was that of an extended 
halo component rotating more slowly than the disk. Here we report the results of H I observations that are 
a factor 5 more sensitive than those of Iswaters. Sancisi. & van der Hulst ( 1997 ). These observations reveal 
that the gaseous halo is much more extended than showed by the previous data. 
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2. Observations 



The present observations were obtained with the Westerbork Synthesis Radio Telescope in the period 
August - December 2002. In total, 20 complete 12-hour observations were performed, using five of the 
standard array configurations. The combination of these different configurations gives a regular sampling 
of the uv plane from the shortest spacing of 36 m to the longest baseline of 2754 m, with an interval of 36 
m. Care was taken not to use data affected by solar interference that might have compromised the detection 
of faint, extended emission from the halo. The effective integration time corresponds to that of 17 complete 
12-hour observations. The observing bandwidth is 10 MHz (corresponding to about 2000 km sec"'), using 
1024 channels (with 2 independent polarizations). An overview of the ob servational parameters is given in 
table |2l The data processing was done using the MIRIAD package (|SaultilTeuben & Wright 19950 . Before 
and after each 12-hr observation a standard calibrator was observed (J2052-I-365 and 3C147) from which the 
spectral response of the telescope was determined. As is standard practice with the WSRT, during each 12- 
hr track no additional (phase) calibrators were observed to monitor the time variation of system properties. 
Instead, the large bandwidth allows to determine these by self-calibration of the continuum image made 
from the line-free channels of the data. This self-calibration was done using a model of the continuum 
emission based on the combination of all observations. An advantage of this approach is that it also gives an 
excellent removal of the continuum sources in the line images, as well as a well-calibrated continuum image 
(see Sec. 3.4). 

Three datacubes were made with three spatial resolutions 23.4" x 16.0", 33.2"x23.9" and 69.6"x58.9" 
(see Table O. The line data were combined and gridded into cubes of 224 channels 8.2 km sec"' wide to 
which additional Hanning smoothing was applied. This results in a velocity resolution of 16.4 km sec '. 
As always with radio observations, the highest-resolution data set has the lowest noise level (0.09 mJy 
beam"'), the 30- and 60-arcsec data sets have noise levels of 0.10 and 0.12 mJy beam"' respectively . The 
3-cr detection limits over one resolution element in the three datacubes are 1.3x10''' cm"2, 6.8x10'^ cm 



and 1.6 X 10 cm . The datacubes were cleaned using the Clark algorithm. In an iterative procedure, 
regions with line emission were identified by smoothing the data to twice the spatial resolution and selecting 
a clip level by eye to define the mask where to clean the data. This procedure was repeated until convergence 
was achieved. 



The H I flux integral is 1.92 • 10 Jy k m sec corresponding to a H I mass of 4.1 • 10^ McQ. The 



mass deriv ed here is 10% more than found bvlSancisi & AUenl (| 19791 ') and by lRupenl (|l99ll '). but the same as 
derived by lBraun. Thilker & Walterbos I (|2003h . 



we assume a distance to NGC 891 of 9.5 Mpc, Jvan der Kruilll98lh ; 1' corresponds to 2.76 kpc 
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3. Observed properties 



3.1. Density map 

In Fig. [T] (right panel) we show the total H I image for NGC 891 at 30" resolution. This image was 
obtained in the standard way by smoothing the datacube to a resolution of 40" and using this smoothed 
cube to create a mask to be applied to the original cube. The mask was produced by blanking emission 
below +3 r.m.s. noise in the smoothed cube. To illustrate the effects of the improvement in sensitivity of 
almost two orders of magnitu de, in Fig. [J we also give the H I density distribution as pub lished in two 
previous studies of NGC 891 (|Sancisi & Allenlll979l : ISwaters. Sancisi. & van der Hulstill997r) . One feature 
immediately clear from Fig. \T\ is that, despite the large improvement in sensitivity, the radial extent of the 
H I disk has not become significantly larger and that in the plane, the new observations do not reveal any 
basically new features (see below). 

However, in the vertical direction the situation is dramatically different and in the new observations 
the H I extends much farther a bove the plane. The extra- planar emission in NGC 891, barely visib le in the 
data of lSancisi & AUenl (119791) and detected up to 5 kpc in lSwaters. Sancisi. & van der Hulsd (119970 . is now 



detected up to a projected distance of more than 10 kpc from the plane everywhere and more than 20 kpc 
in a filamentary structure in the NW. Integrating the flux density located above and below 1 kpc from the 
plane, we find that the HI in the halo represents 29% (1.2 x 10*^ Mq) of the total H I content of NGC 891. 
A s one can see from the lowest contours, t he difference in sensitivity between our total H I map and that 
of ISwaters. Sancisi. & van der HulstI (|1997h is about a factor 5. The H I observations presented here are 
among the deepest ever obtained for an external spiral galaxy. It is therefore quite possible that also other 
galaxies, if observed with comparable high sensitivity, would show similar extended extra-planar emission 
and gaseous halos may be a common feature among spiral galaxies. 

The vertical extent of the H I layer is illustrated in Fig. |2] where we plot the normalised average H I 
column density, based on the 30 arcsecond dataset, for the four quadrants of the galaxy. In the NE, SE and 
SW quadrants the H I halo extends out to about 5' (14 kpc) while in the NW quadrant it extends even further 
out to about 8' (~22 kpc). The density profiles in the NE, SE and SW quadrants are very similar and can be 
modelled quite well with an exponential profile with scaleheight of 50 arcsec (2.2 kpc). In the NW quadrant 
the density distribution seems to follow the same trend out to about 3' (8.3 kpc) after which it becomes 
flatter. This larger extent in the NW quadrant corresponds to a large filament extending out to more than 20 
kpc from the disk (see Fig.[T]). 

Around the NW filament a crowd of high-latitude clouds are observed (Figs. |4] and [T6l ). which are 
probably associated with the filament itself. However, individual gas clouds are also observed at large radii 
and at very anomalous velocities (with large deviations from rotation). The middle panels of Fig. [16] show 
two of these clouds at "apparently" counter-rotating velocities. The first is in the N-W quadrant at velocities 
that differ by about 100 km sec~^ from the velocity of the halo at that position. This cloud is probably 
located in the outskirts of the halo, otherwise the drag force of the halo would invert its motion very quickly. 
The second cloud is detected (Fig.[T6l middle right) at a projected distance of about 28 kpc from the center 
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of the galaxy. These clouds have masses of Mhi ~ 1 - 3 X 10^ Mq. 

As stated above, the new data do not reveal significant new features at large radii in the plane of the 
disk. As the earlier data showed, the H I disk of NGC 891 is not symmetric, being more extended on the 
South-West side. The fall-off and disappearance of the H I disk on the northern side of NGC 891, approx- 
imately coincident with the end of the stellar disk, and the large southern e xtension, confirm the picture of 
lopsid edness already known from previous observations and discussed by iBaldwin. Lynden-Bell & Sancisi 
(Il980l) . However, the fact that, despite the large improvement in sensitivity, NGC 891 does not grow in 
radius but does grow substantially in the vertical direction, poses an interesting question. The presence and 
the origin of outer cut-offs of H I disks have been a matter of debate in the past ye ars. One of the favored 
explanations for their origin has been the effect of the extragalactic radiation field (lMaloneylll993h . In the 
case of NGC 891 such an explanation would immediately encounter a serious objection: why would the H I 
be ionized by a presumably isotropic radiation field in the plane of the galaxy (and account for the northern 
truncation) and apparently not be affected at all in the low-density halo regions? Unless the time scales for 
the disk and for the halo gas and for their respective ionizations are significantly different, it seems more 
likely that the disappearance of H I on the northern side of NGC 89 1 simply marks the outer boundary of 
the gaseous disk. 



3.2. Kinematics 

The kinematics of the H I gas is shown by the position-velocity (p-V) diagrams parallel and perpendic- 
ular to the disk (Figs.[3]and|4ll. In these diagrams, the high-resolution data are shown with thin contours and 
the grayscale, while the thick contours show the low-resolution (60") data at a 3-cr level. The low resolution 
is intended to outline the full extent of the faint emission. The upper plot in Fig. [3] shows the distribution 
along the major axis (i.e. in the plane) of NGC 891. The rotation curve is shown by the sq uares and crosses 



in Fig .m This has been derived by tracing the "envelope', as already done and described by lSancisi & Allen 



(119791) . and using a gas velocity dispersion of 8 km/s. This method is also used by Fraternali (2007, in 
preparation) to derive the rotation curves for the halo gas of NGC 891 at various distances from the plane. 
As already mentioned, the disk of NGC 891 is lopsided, being more radially extended on the receding S-W 
side. Th e rotation in the plane is charac terized by an inner peak produced by a fast-rotating inner ring or 



by a bar (|Garcia-Burillo & Guelinlll995h . This feature is symmetric with respect to the center and also with 
respect to the systemic velocity. Beyond this, the disk is clearly dominated by differential rotation showing a 
roughly flat rotation curve out to a distance of about 6' (~17 kpc) from the center (Fig. [5]). The receding side 
extends further out with an apparent decrease in rotational velocity. However, we do not know the azimuthal 
location of this extension in the plane of NGC 891. If this is not along the line of nodes, the decrease in 
rotational velocity is only apparent and due to projection effects. Note, however, that this extension proba- 
bly represents a large fraction of the outer disk as it can be seen in the p-V diagram of Fig. [5] on the entire 
receding side at projected distances from 0' to 1 1' and on the velocity side closer to systemic. 

Above and below the plane of the disk (lower panels in Fig. [3]), the shape of the p-V diagrams changes 
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dramatically. First, the fast-rotating inner component quickly disappears, indicating that it is confined to 
the inner thin disk. Moreover, the overall shape of the diagram changes from that of a typical differentially 
rotating disk to that of solid-body rotation. In particular, the two diagrams at z = ±3' from the plane clearly 
show the pattern of a slow solid-body rotator. These features are a first indication for a slow rotation of the 
gas above the plane with respect to that in the disk. 

The /j-Vdiagrams in Fig. |4]show the vertical density- velocity structure of NGC 891. The shape of these 
plots is generally tria ngular with the vertex l ocated at the highest rotation velocity (in the plane). This shape 
was already noted by lSancisi & AUenl (119791) in their Figure 7, showing the H I emission on the receding and 
the approaching sides. At increasing distances upward and downward from the plane, the emission tends to 
disappear from the high rotation velocity side and to be restricted closer and closer to the systemic velocity 
(Vsys - 528 km sec~^). This is particularly clear for the emission at very high latitudes. The p-V diagrams 
at distances from R = -1' to R = -3' (i.e. at 1' and 3' from the center of the galaxy on the N-E side) show 
emission at low levels up to about 8' (-22 kpc) and concentrated around the systemic velocity of NGC 891. 
This is the filament visible (Fig. [l) on the NW side of the galaxy. 

The p-V diagrams also show that, if one compares the halo gas in the different quadrants, the kinematics 
of this gas is quite symmetric up to heights of about 3' (~8 kpc). Beyond this, significant differences exist 
between the quadrants. This symmetric kinematics in the lower halo could indicate that the gaseous halo in 
this region is in a equilibrium situation. 

A synoptic view of the kinematical structure of the H I gas in the vertical direction (from the disk 
to the halo region) is shown in Fig. [6] where the thickness of the H I (measured as the first moment of the 
z-distribution) is given as a function of position along the major axis and radial velocity. At the high velocity 
side (close to rotation) the z-distribution is much thinner than at the lower rotational velocity side as also 
seen in Fig. |4l 



3.3. UGC 1807 

The present observations show also UGC 1807, a small, gas-rich companion located at a projected 
distance of about 30' (~ 80 kpc) from NGC 891 (Fig. IT]) (not discussed in previous papers on NGC 891). 
In the optical, UGC 1807 appears as a LSB galaxy oriented almost face-on. Our H I data show a regular 
distribution of H I with a symmetric velocity field (Fig. IT]). The outer radius of the H I disk is 4 kpc and the 
total H I mass corrected for the primary beam attenuation is 4.5 x 10^ M©. We note that at the position of 
UGC 1 807 this correction is quite large, about a factor 6. 

The kinematic parameters of UGC 1 807 have been derived with a tilted-ring fit of the velocity field. 
We found a systemic velocity of 627 km sec"', 100 km sec"' larger than that of NGC 891, and an inclination 
of about 15°± 5°. The rotation curve (Fig. IT) is very regular and shows a slow rising in the inner regions 
up to a rotation velocity of about 100 km sec"' at the last measured point. However, because of the low 
inclination and its large uncertainty, the amplitude of the rotation curve is very uncertain: it can be as low 
as 70 km sec"' or as high as 140 km sec"'. Using 100 km sec"' and the measured outer radius, we find for 
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UGC 1807 a total mass of 9 x lO'^ M, 



3.4. Radio Continuum 



Thanks to the broad band and the long integration time of the present observations we have been 
able to construct a very deep radio continuum image (F ig. [8]l, similar in quality to the 21-cm continuum 
image published by lDahlem. Dettmar & Hummell (119941) . The spatial resolution is 17 x 12 arcsec, the RMS 
noise is 23 juJy beam"^ Our data do not contain, however, the information on the polarization of the radio 
emission. In Fig.|9]we show the normalized ver tical brightness distribution of t he rad io continuum halo in the 
four quadrants of NGC 891. As was found by lDahlem. Dettmar & Hummell (119941) . the vertical brightness 
distribution of the radio continuum closely follows an exponential profile with a scaleheight of 25 arcsec 
(corresponding to 1.15 kpc) in the NE, SE and SW quadrants and a bit larger (1.3 kpc) in the NW quadrant. 
In the SW quadrant the radio continuum halo seems to be somewhat more extended. As in the earlier data on 
NGC 891, there is also a strong N-S asymmetry of the continuum emission in the disk. Dahlem et al. argue, 
based on the similarity in the distribution of the Ha emission and the radio continuum, that both are the 
result of outflows from the disk driven by star formation. In Fig.|9]we also indicate the vertical distribution 
of the Hof emission of the halo, as determined by Rand & Heald (priv. comm.). 



4. Models 

Here we investigate the structure and kinematics of the extra-planar H I emission in NGC 891 by 
constructing model datacubes and corn paring these with the observations. We first consider, following 



Swaters. Sancisi. & van der HulstI (|l997l ). basic models (Section |4~T]) in which the extra-planar emission is 
produced by a strong warp, a flare and a co-rotating thick disk. All these models are easily ruled out by 
the comparison with the data. Then, in Section 14.21 we consider more sophisticated models made of two 
components (disk -i- halo) with different kinematics. 



4.1. Basic models 



We build model datacubes using a modified version of the GIPSY (|van der Hulst et al.ll 19921 ) program 
GALMOD. This program assumes axi-symmetry and an H I layer made of concentric rings. For each ring 
we define the H I column density, the thickness of the layer and the geometrical and kinematic parameters. 
The H I disk of NGC 891 is not symmetric, being more extended on the SW (receding) side of the galaxy. 
In order to avoid the complications arising from this lopsidedness we construct the models only for the NE 
(approaching) side of the galaxy. 

We start with a thin disk model. The H I radial density distribution in the plane of the disk is derived 
by considering only the NE side of the galaxy and the region within \z\ < 30". The result is plotted in Fig. 
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[TOl We approximate the observed distribution with an exponential law at large radii and a depression in the 
inner regions as described by: 



Z(/?)-Eo^l + — j exp (-/?//?.), (1) 

where Sq is the central surface density, a is an exponent defining the tapering towards the center and /?* is 
a scale radius such that the peak of the distribution is located at R = R^ia - 1). The parameters of the fit 
shown in Fig.[lO]are Eq = 6.2 x 10~^ Mopc"^, a = 7.8 and /?* = 1-2 kpc. In the inner regions of the galaxy 
(/? < 3 kpc) this model fails to reproduce the inner ring that is clearly visible in Fig. [5] To reproduce also 
this inner ring, we add an extra exponential component to our model. The parameters of this component are 
So = 6.3 MqPC~^ and R^ = 1.2 kpc. The approaching and receding rotation curves in the disk of NGC 891 
are shown in Fig. [5] In our models, we use the curve derived for the N-E (approaching) side of the galaxy 
and we assume a velocity dispersion o"gas = 8 km sec"^. 

NGC 891 has a mild warping of the outer disk, as can be seen in the top-right plot in Fig. [TT] In all 
models, we have included this mild warp by varying the inclination and position angles of the rings beyond 
R = 13 kpc. The warping angles that fit the data are quite small (1.5° in p. a. and 4° in inclination). The 
inclination angle of the inner disk is taken to be 90°. The vertical distribution has been modeled with an 
exponential profile with scaleheight h^isk - 0.2 kpc. 

Fig.[lT](leftmost column) shows the best fit thin disk model for NGC 891. The rightmost column shows 
the data channel maps at full resolution. An important property of the data is that the apparent vertical 
thickness increases from very thin in the Vhei = 292 km sec"' channel (corresponding to the extreme rotation 
velocity) to a broad vertical distribution in the bottom channel (close to the systemic velocity, Vgys = 528 
km sec"'). This is a major constraint for all our models as it indicates that there is no halo gas with rotation 
velocities as high as those of the gas in the disk. Clearly the thin disk model only reproduces the top channel 
at Vhei = 292 km sec"' . The gas visible in this channel is likely to be at the line of nodes and its thickness 
is a measure of the intrinsic thickness of the thin disk. However, because of the insufficient resolution 
(HPBW = 19" =^ 0.9 kpc in the vertical direction) this channel gives only an upper limit to the thickness of 
the thin disk (see below). 

The second column of Fig.[TT]shows the effect of a strong warp along the line of sight in addition to the 
mild warp described above. Such a strong warp model has been propose d by BecQuaert & Combes (|l997|) 
to explain the extra-planar emission observ ed in NGC 891 with the VLA ( Runenlll991 ), in contrast with the 
model of a thick, lagging layer favored by ISwaters. Sancisi. & van der Hulst "(Il997h . In the present model 
we have introduced a linear change in the inclination of the outer rings (beyond R = 13 kpc) that reaches 
a tilt of 25° in the outermost ring with respect to the inner parts (i = 90°). This model comes close for the 
top and the bottom channels of Fig. [Ubut totally fails to reproduce the middle ones. The characteristic 
feature of any line-of-sight wa rp model is the "bu tterfly" opening visible in the middle channels (around 
Vhei = 366.2 km sec"') (see also Gentile et al.ll2003 ). Such a pattern is totally absent in the present data and, 
therefore, we can rule out that the extra-planar emission in NGC 891 is produced by a line-of-sight warp. 

The middle column of Fig. [TT] shows the effect of a large flaring of the H I layer. In this model, the 
thickness of the H I layer in the inner regions is as small as in the thin-disk model (/Jdisk - 0.2 kpc) but 
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it increases outwards and reaches /jjisk - 3 kpc in the outer rings. Such an unrealistically large value is 
necessary to produce the emission observed at high latitudes. This model generates some channel maps 
fairly similar to those observed (e.g. at Vhei = 366.2 km sec"'), but it completely fails to reproduce the thin 
structure in the top channel maps. This is because in the model the gas in the outer flare rotates as fast as the 
gas in the plane. Also this model, therefore, can be ruled out. As already mentioned above, from the thin 
structure observed at the highest rotational velocities (top Fig. [TTI) it is possible to derive an upper limit to 
the thickness of the thin disk and to set an upper limit to its flaring. We find that the inner scaleheight of the 
disk is less than /jdisk < 0.3 kpc, while in its outer regions the disk can flare up to at most husk ~ 0.5 kpc. 

We now consider models made of two components: 1) a thin disk (like that in the leftmost column in 
Fig. [TTI ) and 2) a thicker layer or halo. For the vertical gas distribution of the latter we use a density profile 
described by the following empirical function: 

sinh(z//!haio) 

COSh(z//jhalo)^ 

where z is the vertical coordinate, is the surface density in the plane and /ihaio is the halo scaleheight. With 
this formulation, the vertical density of the thick layer is zero in the plane, then rises reaching its maximum 
at z = 0.88/thaio and declines nearly exponentially further out. The half width at half maximum (HWHM) of 
this distribution (defined as the z where ^(z) = ^(0.88/zhaio/2) is at z = 2/jhaio- Between disk and halo there is 
a gradual transition. In this way the spatial coexistence between the two components is minimized: the halo 
takes over when the disk is fading out. 

In order to model the gas density in the halo, we have extracted and de-projected the radial distributions 
at various heights. The left panel of Fig. [12] shows (squares) two of these deprojected radial distributions 
atz = 2.8 and z = 5.6 kpc (average of the NE and NW quadrants). The shapes of these distributions change 
with the height from the plane: as the distance from the plane increases, the radial distribution becomes 
flatter and less concentrated to the center (cf. Fig.[TOl). We have fitted the data with the function: 

\ v^m^*^^''^halo(^)) 
Phalo(^,Z) ^ ^{R) ; (3) 

where S(/?) and ^(z; /ihaio) are given by equations 1 and 2 but the scaleheight /zhaio varies with R. Fig. [12] 
(right) shows the fitted values (squares) of /jhaio as a function of R and a power law fit (solid line). With 
this parametrization the scaleheight of the halo varies from /jhaio - 1-25 kpc (HWHM= 2.5 kpc) in the 
central regions to /ihaio ~ 2.5 kpc (HWHM~ 5 kpc) in the outer parts. The other parameters of the fit are: 
Xo = 1-4 X 10"' Mopc"^, a - 4.5, /?* - 1.9 kpc. The final result is shown (hnes) in the left panel of Fig. [12] 
(at z = 2.8 and z = 5.6 kpc from the plane). 

The forth column in Fig. [TT] shows a two-component model in which the halo co-rotates with the disk. 
The mass of the halo component is Mhaio = 1-25 x 10^ Mq, 30% of the total H I mass. Clearly, this model 
does not correctly reproduce the channel maps at high rotation velocities (top two panels) which appear 
much thicker in the model than in the data. 
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4.2. Lagging-halo models 

The failure of the two-component model with a co-rotating halo to reproduce the channels near the 
extreme rotation velocity suggests that the gas above the plane is rotating more slowly than that in the 
pla ne, i.e. the halo is lagging in rotation w ith respect to the disk. Evidence for su ch a lag was found 



bv iSwaters. Sancisi. & van der Hulsd (119971) and is also observed in the ionized gas (iHeald et all I2006bl : 
Kamphuis et al.ll2007n . Fig. [13] shows four models with lagging halos where the rotation velocity decreases 
with increasing height above the disk. All these models consist of a thin disk (see Fig. [TT]) and a halo 
component with the density distribution described above. Table |4] lists the kinematic parameters used in the 
lagging-halo models. The first model (leftmost column) is that of a lagging halo with the vertical gradient in 
rotation velocity independent of radius. The rotation curve in the disk is that shown in Fig. [5] (neglecting the 
inner fast-rotating ring). We have adopted, after a few trials, a constant negative vertical gradient in rotation 
velocity of Avrot/Az = -0.55 km sec^^ arcsec"^ - -12 km sec"' kpc"' for the halo component. 

It is clear that this simple lagging-halo model reproduces the main features of the data much better than 
the previous models. The structure in the upper channel maps is as thin as in the data and it becomes thicker 
as one approaches the systemic velocity (bottom row). 

However, this model is still not completely satisfactory. In particular, near the systemic velocity, the 
radial extent of the halo gas is much narrower than in the data. One way to improve the model is to increase 
the velocity dispersion of the halo gas (second column of Fig.[T3]). For this model all values of the parameters 
have been kept the same as in the previous model, except the velocity dispersion of the halo gas which has 
been increased to cThaio = 25 km sec"' and the halo rotation which has been decreased. The higher velocity 
dispersion is physically plausible since the gas in the halo is expected to be kinematically "hotter" and to 
have a more complex motion than the gas in the plane. However, to such increasing velocity dispersion 
corresponds a decreasing halo rotation and, therefore, an increasing vertical velocity gradient to Avrot/Az = 
-0.8 km sec"' arcsec"' - -17.4 km sec"' kpc"'. Fig. [T3] shows that increasing the velocity dispersion of 
the halo gas indeed improves the match of the radial extent of the halo H I in the channel maps close to the 
systemic velocity (bottom row). 

An alternative approach to the increase of the velocity dispersion is to introduce a systematic non- 
circular motion in addition to the rotation of the halo gas. Such non-circular motions have been found in 
similar studies of other galaxies. Large-scale inflows toward the galaxy center have been dis covered in 
NGC2403 (vrad - -15 km sec"') JFratemaU et allboOll) and in NGC 4559 JBarbieri et al.ll2005l) . The third 



column of Fig. [T3] shows the eff'ect of an overall radial motion of the halo of |vradl ~ 25 km sec"'. This is 
of about the same amplitude as found in NGC 2403 (although in an edge-on galaxy one cannot discriminate 
between in- and outflow). The models also show that the effects of a higher velocity dispersion or radial 
motions in the halo are very similar. It is, therefore, not possible to discriminate between an in/outflow and 
a higher velocity dispersion of the halo gas. 

All the models described above do reproduce most of the features present in the data, but do not fully 
account for the shape of the middle channel maps of NGC 891. Let us focus in particular on the channel map 
at Vhei = 366.2 km sec"'. The shape in the data channel map is roughly triangular, while in the models it is 
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more boxy. Clearly, in the model channels near Vhei - 366.2 km sec" , there is too much halo emission near 
the center. In order to obtain the triangular shape, we consider two possibilities. The first is that the inner 
regions of the halo are depleted of gas (much more than shown in Fig. [TOl ). This obviously would decrease 
the amount of gas near the center in the models. The second possibility is that the gradient in rotation 
velocity in the inner parts of the galaxy is larger than in the outer parts. First we consider the possibility 
of a stronger central gas depletion in the halo. Since the gas density in the halo was derived from the data 
(without any assumptions about the kinematics), and modeled accordingly, we do not expect the used H I 
density to be significantly different from the actual one. The main source of errors is the de-projection of 
the radial profiles. However, this is not expected to have a strong effect beyond R > 3 kpc, while in order 
to reproduce the channel map at Vhei = 366.2 km sec"' the gas density should be significantly lower than 
the one used here for radii as large as /? ~ 8 kpc. Therefore, it seems unlikely that a central depletion is the 
explanation. 

Consider instead the second possibility. The forth column from left in Fig. [13] shows a model (named 
shallow rise) in which the vertical gradient in the rotation velocity varies with R. In the inner regions, 
the gradient is quite large (Avj-ot/Az = -43 km sec"' kpc"' at R = 0), and it decreases linearly with R 
(A(Avi.ot/Az) /AR ~ 2.5 km sec"' kpc"^) until it reaches a value of Aviot/Az = -14 km sec"' kpc"' and 
remains constant further out. The shape of the rotation curve of the halo changes with distance from the 
plane: its inner rising part becomes shallower and shallower. This behavior is illustrated in Fig. [15] Clearly, 
this model gives a better representation of the channel at Vhei = 366.2 km sec"' while keeping almost 
unchanged the other channel maps (note that this shallow-rise model has also a high cr for the halo gas and 
an inflow similar to the previous two models). The rotation curves for the halo gas at various distances from 
the plane, derived directly from the observations, will be presented and discussed by Fratemali (2007, in 
preparation). 

We further compare the models by inspecting the position-velocity cuts perpendicular to the major axis 
of the galaxy. Fig. [T4l shows two sets of such cuts taken on the NE side of the galaxy at a distance of 1' (2.8 
kpc, bottom) and 2.7' (7.5 kpc, top) from the center. The emission shows a characteristic triangular shape. 
The thin-disk is visible at all velocities between systemic (vhei - 528 km sec"') and the maximum rotation 
(vhei ~ 290 km sec"'), whereas the halo appears to have its maximal extent near the systemic velocity. It is 
clear that the shallow-rise model best reproduces the triangular shape as well as the difference between the 1' 
and 2.7' plots. The white squares in the data plots show the rotational velocities adopted for the shallow-rise 
model. 

The comparison of the various models leads to the conclusions that the thin disk of NGC 891 is sur- 
rounded by an extended gaseous halo which is rotating more slowly than the disk and contains almost 30% of 
the neutral gas of the galaxy. The kinematics of the gas in the halo can be best explained by assuming a verti- 
cal gradient in the rotation velocity. Moreover, this gradient is stronger in the inner regions (Avrot/Az ^ -43 
km sec"' kpc"') than at larger radii (Avrot/Az ^ -14 km sec"' kpc"'. Finally, the velocity dispersion of the 
halo gas is higher than that of the gas in the disk (crhaio - 20 - 25 km sec"' vs CTdifik - 8 km sec"') and/or 
there are significant radial motions in the halo gas (|vradl ~ 25 km sec"'). 
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5. Discussion 

In the previous Section we have studied with 3-D models the extended extra-planar H I emission of 
NGC 891 and concluded that this galaxy has a massive halo of neutral gas rotating more slowly than the 
disk. Here we discuss the properties of this H I halo, its possible origins and the comparison with data at 
different wavelengths. 



5.1. The structure of the H I halo 



The H I halo of N GC 891 is the most extended and massive of those found to date in a spiral galaxy 
(iFratemali et aLll2007al) . The H I mass above z = I kpc is about 30% of the total H I mass. The distribution 
of gas in the halo appears fairly symmetric and regular in the four quadrants of the galaxy up to about 8 kpc 
from the plane (Figs. [J and O. In three quadrants, the H I extends further up to about 14 kpc, whereas the 
N-W quadrant is dominated by an extended filament reaching up to 8' (~22 kpc). Radially, the halo extends 
to the end of the disk on the N-E side but stops earlier on the S-W side where the disk is more extended. 
This may be an indication that the halo is closely connected to the inner disk of NGC 891. 

Also the kinematics of the halo is symmetric and regular in the four quadrants up to a height of about 
8 kpc (Fig. |3]l. The gradient in rotation velocity can be measured up to about 5 kpc. Above that the rotation 
velocity continues to decrease further, but the gas density is too low to derive reliable rotation velocities 
(IFratemah et al.ll2005r) . The gas in the halo may have radial motions and/or higher velocity dispersion than 
the disk gas, as the above model ana lysis indicates (Section 4.2). Vertical motions may also be present (as 
in NGC 2403, IFratemali et al.l (|200lh ) but because of the inclination of the galaxy they cannot be observed. 



The halo of NGC 891 shows individual features (streams and compact clouds) somewhat "separated" 
either in location or kinematics from its main bulk. The most prominent feature is the extended filament in 
the N-W quadrant (see Figs. [T] and |4ll. This feature has a mass of more than 1.6 x 10^ Mq, it extends out to a 
projected radius of about 10 kpc from the center of the galaxy and vertically up to 22 kpc (Figs. l4l[T6l). The 
full length of the filament is about 30 kpc, its velocity width is Ay ^ 100 km sec"' and its mean velocity is 
close to systemic for z > 10 kpc. 

A key question is whether the H I halo of NGC 891 is a diffuse medium (filling factor ~ 1) or it is en- 
tirely made of discrete compact clouds. A first inspection of the datacube seems to indicate that most of the 
gas belongs to a smooth and coherent, differentially rotating structure. This picture of a "diffuse medium" 
may, however, be erroneous and be produced, at least to some extent, by projection effects along the line of 
sight. Only clouds with anomalous velocities or located at large distances from the plane would stand out 
clearly (such as those shown in Fig.[T6l). These clouds, together with the filament on the NW quadrant may 
be evidence that the halo of NGC 891 is indeed made, at least partly, of individual gas complexes. These 
complexes would be similar to those observed in non-edge - on galaxies like NGC 2403 and NGC 69 46 where 

More- 



projection effects are less important (|Fratemali et al ] l2002l : Isoomsma et al.ll2005al : 



Boomsma 
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over, they would have similar properties as the HVCs and IVCs of the Milky Way (jwakker & van Woerden 
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19971) and the clouds seen near M31 dWestmeier. Braun. & Thilkeill2005h . A structure like the long NW 



filament would probably appear as a "Complex A" or "Complex C" of the HVCs to an observer inside 
NGC 891. 



5.2. Origin of the H I halo 

What is the origin of the H I halo of NGC 891? We discuss here two possibilities: the galactic fountain 
and the accretion from outside. The main observational facts to be accounted for are: 1) the halo is very 
extended and massive (about 30% of the total H I mass); 2) the halo kinematics is dominated by differential 
rotation and the rotation velocity decreases with height from the plane; 3) a high velocity dispersion or non- 
circular motions (maybe inflow?) are also present; 4) the halo is structured in clouds and filaments, some of 
which are at very anomalous (counter-rotating) velocities; 5) As the distance from the plane increases, the 
radial distribution of the gas in the halo tends to become flatter and less concentrated to the center than in 
the disk. On the southern side, the halo is radially less extended than the disk. 



5.2.1. Fountain 



There is considerable evidence from Ha radio continuum and X-ray data pointing at a galactic fountain 
mechanism playing a majo r role in forming the gaseous halo of NGC 891. The Ha image of NGC 891 
(|Dettmanll990l : lRandlll990l ) indicates that the star formation rate in the disk, especially on its northern side, 
is very high. There is strong radio continuum emission (thermal and non-thermal) in the disk and also (non- 
thermal) in the halo, extending up to ~10 kpc from the plane (see Fig. [H). This extent is close to that of the 
H I halo and indicates the presence of magnetic fields and relativistic electrons in corresp ondence of the H I 
Eas^ 



Also, there is a clear correlation between this radio emission and the Ha halo (|Dettmarlll990l : 



Rani 



1990l ; lDahlem. Dettmar & Hummelll 19941) . In particular, the northern side of the galaxy is much brighter in 
both components and in the NW quadrant the radio halo has a larger scaleheight (Fig. |9l). A corresponding 
N-S asymmetry is also seen in the H I halo. Fig.[6l which shows the thickness of the H I layer as a function of 
position and velocity, indicates that in the inner regions of the disk (corresponding to small R and velocities 
furthest away from systemic), the H I disk is thicker on the northern (i.e. approaching) than on the southern 
side. The northern side of the disk is where star formation and consequently the fountain are strongest. 

The ionized gas in the halo of NGC 891 h as a smooth component (Diffuse Ionized Gas) but it also shows 
filamentary structures iHowk & Saya ge^^200oV CoUimated Ha filaments are seen reaching up to more than 
2 kpc from the plane (IRossa et all 12004) . Moreover, the radio emission is strongly polarized. Both these 
features suggest the presence of a uni form magnetic field in the halo, ind icating that the outflowing gas has 
also torn the B-field out of the disk (IDahlem. Dettmar & Hummelll 19941) . Interestingly, also the long H I 
filament found in our data is similarly oriented. However, in the Ha images there is no trace of an ionized 
counterpart of the H I filament to the detection limit. 

Hot coronal gas was first revealed in NGC 891 by ROSAT jBregman & Pildisll 19941) . Recently, Chandra 
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and XMM data have confirmed this detection and have shown the pr esence of filamentary substructures in 
the halo component extending up to about 5-6 kpc from the plane dStrickland et al.ll2004l) while the soft 
halo component seems to be concentrated to the inner disk (ITuUmann et alj|2006al) . These authors were 
able to show for a sample of galaxies, including NGC 891, that the amount of hot gas is proportional to 
the mechanical feedback from supemovae, concluding that this hot gas is alm ost certainly produ ced by a 
fountain. The X-ray emission also correlates fairly well with the Ha emission (|Rossa et alJl2004l) . and the 
X-ray spectrum of the halo emission in the range 0.3-2 keV is fitted by a plasma with a temperature of 0.23 
keV (3.7 X 10^ K). 

Besides the H I there are in NGC 891 also other "cold" extra-planar components suc h as dust and CO. 
The dust sho ws up as absorption features against the star light up to heights of about 2 kpc (IHowk & Savage 



19971 . l2000h . Such features have masses of about 10^ Mq (using the galactic dust-to-gas ratio). This 



cold medium d oes not seem to correlate, on small spatial scales, with the warm medium observed in Ha 
(|Howk & Sava ge 2000; Rossa et al. 2004); thus the two phases ap pear to be physically distinct. CO emission 
is als o observed up to about 2 kpc from the plane of NGC 891 ( Garcia-Burillo et al 1992 : Sofue & Nakail 
1 1993k This extra-planar component of molecular gas comprises about a 35% of the total amount of molec- 
ular gas. 

All this evidence suggests that a classical galactic fountain (ionized outflow, cooling, and cold inflow) 
is active in NGC 891 and is responsible for much of the gaseous components, at least in the lower halo (z ^ 5 
kpc). These components have a total mass of the order of a few x 10^ Mq. Among them, the H I halo is 
the most massive and the most extended reaching a distance from the plane of more than 20 kpc. This latter 
may be the result of an observational bias or the indication that a diff"erent mechanism, perhaps accretion 
from outside (see below), is responsible for the formation of the upper layers. 

Galactic fountain models can easily explain the distribution of the HI in the halo of NGC 891, with the 
exception of the exte nded filament (see below). Indeed, only a few percent of the energy from supernovae 
would be sufficient (|Fratemali & Binneyll2006h . However, such models fail to reproduce the kinematics 
of the halo g as (both ioni zed and H I). The ionized component has been initiall y studied using long -slit 
spectroscopy (|Randlll997n and more recently using the WIYN integral field unit (|Heald et al.ll2006br) and 
TAURUS spectroscopy (|Kamphuis et al ] l2007h . The gradient in rotation velocity derived for the ionized gas 
in the latter study is in good agreement with that observed in H I about -15 km sec~'kpc~' for 1 < z < 5 
kpc. This result indicates that the H I and Ha phases are part, at least in the lower halo, of the same 
phenomenon and, therefore, probably also have the same origin. A similar coupling between ionized an d 
neutral halo gas has also been found in the spiral galaxy NGC 2403 (IFratemali. Oosterloo & Sancisill2004r) . 
This ob served gradient in rotation velocity with z is significantly higher than predicted by the fountain 
models (|Fratemali & Binneyll2006l : iHeald et al.ll2006bh . The only way to reconcile these models with the 
data is to find a mechanism that would make the fountain gas loose part of its angular momentum. One 
possibility is that it is braked by the hot halo. However, in this c ase, the angular momentum transfer to the 
hot halo would speed it up very quickly (|Fratemali et al.ll2007br) . A second, more promising possibility is 
that the fountain gas interacts with low angular momentum material accreted from the IGM (Fratemali & 
Binney, in preparation). 
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5.2.2. External origin 

It is clear from the foregoing discussion that a "fountain" origin for most of the H I in the halo of NGC 
891 is very likely. But a continuous supply of low angular momentum material from outside may also be 
needed to account for the observed H I kinematics. The amount of accreted gas, however, can be small, 
perhaps not more than 10 percent of the gas present in the halo (Fratemali & Binney, in preparation). 

Is there any direct observational evidence of such gas infall in NGC 891? An obvious candidate for 
accretion is the filament in the NW quadrant extending up to 22 kpc from the plane. As mentioned above, 
this filament has a projected velocity close to systemic. This may indicate either of these two possibilities: 1) 
the filament belongs to the inner parts of the halo and has a very small z component of the angular momentum 
or 2) it is located in the far outer parts of the halo. For the filament to be produced by a galactic fountain, one 
would re quire very high kick velocit ies and a very high initial kinetic energy. Using the galactic fountain 



model of iFratemali & Binneyl (120061) . we estimate that the required kick velocities to send material up to 
z = 15 kpc range from 240 km sec~^ for a starting radius of 10 kpc to 425 km sec"' for a starting radius 
of /? = 2 kpc. For a mass of the filament of 1.6x10' Mq, the kinetic energies would be 1 - 3 x 10^^ 
erg. Therefore, to produce such a filament in one single event, hundreds of thousands of supernovae would 
be required. This seems quite unlikely. If instead the gas complex had an external origin, its filamentary 
structure and its kinematics might be more easy to explain: a cloud falling onto NGC 891 would be stretched 
out and might speed up to full rotation as it gets close to the disk. 

An external origin for the filament may be as cold accretion from the IGM. Or it may be through a 
condensation of the hot coronal gas, possibly triggered by the galactic fountain. Alternatively, the filament 
could be the result of the merging of a companion with NGC 891. The stellar part of the companion could 
be elsewhere or too faint to be seen and the filament would be the only prominent remnant of the process. In 
this respect we note that deep H I observations are a powerful tool to look for this kind of events and could 
be used to detect merging at much larger distances than optical data. 

Finally, we note that the filament roughly points at the projected position of the companion UGC 1807 
(Fig. IT]) and one may wonder whether it could have originated from a close encounter of UGC 1 807 with 
NGC 891. The filament would have to be either gas pulled up from the outer disk of NGC 891 or gas 
stripped from the companion. This is not impossible, in spite of the present large separation of UGC 1807 
from NGC 89 1 and of its symmetrical structure and kinematics, which would seem to exclude recent tidal 
encounters. The projected distance of about 80 kpc could have been covered by the companion with a 
velocity of 100-200 km sec"' in 8-4 x 10^ yrs. This is sufficiently longer than the dynamical time scales 
in the outer parts of the companion (2 x 10^ from 100 km sec"' rotation at 4 kpc radius) to account for the 
absence of any sign of disturbances. On the other hand, the filament, presumably located in the outer parts 
of NGC 891, would have survived long enough to be observed now. In this connection it is interesting to 
draw attention to another puzzling, long-known structure revealed by the H I observations of NGC 891: the 
extended southern tail. This is difficult to maintain for a long time against the efi'ect of differential rotation 



(|Baldwin. Lynden-Bell & Sancisilll980h and the possibility of a recent origin (less than 1 x 10^ yrs) would 



therefore be quite attractive. It is therefore interesting to ask whether a tidal encounter with UGC 1807, as 
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considered here for the origin of the filament, could also be responsible for the north-south H I lopsidedness 
of NGC 891. Note that the mass of the companion (-10% of that of NGC 891) is not negligibly small. 

Other indications of accretion are the counter-rotating clouds shown in Fig. [16] These clouds cannot 
be produced by a galactic fountain as they have an orientation of the angular momentum opposite to the gas 
in the disk. They are almost certainly located in the outer halo (one is observed at a projected distance of 28 
kpc), otherwise the drag by the co-rotating H I would have quickly inverted their motion. Note that one of 
these cloud complexes is spatially (and kinematically) very close to the filament. The total H I mass present 
in those accreting clouds (filament included) is of a few 10' Mq. Considering a typical inf ailing time of the 
order of 10^ yr, this would correspond to an accretion rate of « 0. 1 Moyr~' . This latter should be considered 
as a lower limit to the total accretion rate as most of the accreting material may be impossible to detect after 
it has interacted with the fountain halo. 



6. Summary of results 

The deep H I observations of NGC 891 reported here have shown the following: 

1. The H I halo is considerably more extended than shown in previous observations. It is detected 
out to a vertical distance of about 22 kpc from the galactic plane and contains 1.2 x 10^ Mq amounting to 
about 30% of the total H I. The halo is not smooth. It shows structures on various scales and there are, in 
particular, also H I clouds with very anomalous (counter-rotating) velocities. These have masses of about 
1 X 10^ Mq. Furthermore, there is a striking long filament extending (in projection) up to more than 20 kpc 
from the disk and containing more than 1.6 x 10' Mq. 

2. The H I halo has overall regular differential rota tion but it rotates more slowly than the disk (as 



already found by lSwaters. Sancisi. & van der Hulsli (|l997h '). The vertical gradient of the rotational velocity 



is about -15 km sec~^ kpc~^. The shape of the halo rotation curves does not remain constant with distance 
from the plane, but its rising part becomes shallower and shallower with height. Random motions in the 
halo or systematic deviations from circular motion (in- or out-flows) may also be present. 

3. A significant fraction of the H I halo must come from a galactic fountain, as the high star formation 
rate in the disk strongly suggests. The radio continuum, the Ha and the X-ray data are all corroborating 
evidence. However, the kinematics of the halo suggests the need of an interaction between the fountain gas 
and low angular momentum material. Such material may be supplied by gas accretion from the surrounding 
IGM. The long H I filament and the counter-rotating clouds may be direct evidence of such accretion. 
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Table 1. Optical and radio parameters for NGC 891. 





NGC 891 


ref 


Morphological type 


Sb/SBb 


1,2 


Radio continuum centre {a, S J2000) 


2'" 22" 33.0" 


5 




42° 20' 57.2" 


5 


Distance (Mpc) 


9.5 


1 


Lb (Les) 


2.6x10"^ 


3 


Scale length stellar disk(kpc) 


4.4 


4 


SFR(Moyr-') 


3.8 


6 


Systemic velocity (km s" ' ) 


528 ± 2 


5 


Total H I mass (Mq) 


4.1 X 10' 


5 


H I inclination (°) 


> 89 


5 


Mean P. A. (no warp) (°) 


23° 


5 


Total mass (Mq) 


1.4 X 10" 


5 



References. — (1) Ivan der Kr ui? h98l b, (2) 

iGM-c ia-Burill o & Guelii] 1 19951), ( 3) |de Vaucoule urs et al J 
TT99lh. (4) IShaw & Gilmord h989h . (5) This work, (6) 
|PopescuetalJj2004l) 



Table 2. Observational parameters for NGC 891. 



NGC891 



Observation dates Aug. - Dec. 2002 

Total length of observation (hr) 20 x 12 

Effective integration time (hr) 204 

Pointing R.A. (J2000) 2^ 22" 33.76" 

Pointing Dec (J2()(X)) 42° 20' 57.1" 

Velocity centre of band (km sec"') 800 

Total bandwidth (MHz) 10 

Total bandwidth (km sec"') -2000 

Number of channels in observation 1024 

Channel separation in observation (kHz) 9.77 

Number of channels of final datacube 224 

Channel separation of final datacube (kHz) 38.5 

Velocity resolution (after Hanning) (km sec"') 16.4 
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Table 3. Parameters of the data cubes. 





FuU resolution 


Intermediate resolution 


Low Resolution 


Spatial resolution (") 


23.4 X 16.0 


33.2 X 23.9 


69.6 X 58.9 


P.A. of synthesized beam (°) 


12.3 


16.0 


16.0 


Beam size (kpc) 


1.08 X 0.74 


1.53 X 1.09 


3.2 X 2.7 


R.m.s. noise per channel (mJy beam"') 


0.090 


0.10 


0.12 


Minimum detectable column density (3o"; cm"^) 


13.0 X 10'* 


6.8 X 10"* 


1.6 X 10'** 


Minimum detectable mass (3cr; Mq per resolution element) 


9.3 X lO"* 


1.0 X 10' 


1.3 X 10^ 
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Table 4. Parameters of the lagging -halo models shown in Figs. [T3] and [141 



Model 


AVrot/Az 

(km sec"' kpc"') 


O-gas 

(km sec"') 


Vrad 

(km sec" ' ) 


Lagging halo 


-12.0 


8 





High o" 


-17.4 


25 





Inflow 


-13.0 


8 


-25 


Shallow rise 


-43.0 ^ -14.0 


20 


-15 
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Fig. 1. — Left: total H I image as obtained from the observations described in this paper. The outer grey 
contour comes from the 60" data and its level is 0.005 • 10^' cm~^. The black contours come from the 30" 
data with contour levels 0.01, 0.02, 0.05, 0.1 (thick hne), 0.2, 0.5, 1.0, 2.0 and 5.0 • 10^' c m~^. For compar- 
ison o f the total H I images of NGC 891 published over the years, we show the one from lSancisi & Allen 
(119791) (top right) with spatial resolution of about 30" and contou r levels of 0.5, 1.0, 1.5, 2.4, 4.1, 5 . 7, 7.3 , 
8.9 and 10.5 • 10^' cm"^ and the total H I image as published by ISwaters. Sancisi. & van der Hulsd (Il997h 
(bottom right). The spatial resolution of the last image is about 20" and the contour levels are 0.07, 0.17, 
0.46, 1.1, 2.3, 4.1, 6.4 and 9.2 ■ 10^' cm~^. To aid the comparison, two stars labeled 'A' and 'B' are marked 
in each figure 
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Hl z distribution 




z (arcmin) 

Fig. 2. — Normalized average H I column densities along the z direction for the four quadrants of NGC 891. 
The spatial resolution of the observations is indicated. (1' corresponds to 2.76 kpc) 
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Distancc along major axis (arcmin) 




Distance parallel to major axis (arcumi) Distance parallel to major axis (aiemin) 



Fig. 3. — Position- velocity plots parallel to the plane for NGC 891 made from the high-resolution data cube 
(thin lines and greyscale) and the 60-arcsec data (thick hne) Contour levels are in steps of l.Scr for the 
high-resolution data and at 3cr for the low-resolution data. 
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Helioqentnc velocity (km/s) Heliocentric velocity (km/s) 



Fig. 4. — Position- velocity plots perpendicular to the plane for NGC 891 made from the high-resolution data 
cube (thin lines and greyscale) and the 60-arcsec data (thick line). Contour levels are in steps of 1.5cr for 
the high-resolution data and at 3cr for the low-resolution data. 
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Fig. 4 (cont'd). — Position- velocity plots perpendicular to the plane for NGC891 made from the high- 
resolution data cube (thin lines and greyscale) and the 60-arcsec data (thick line). Contour levels are in steps 
of 1.5o" for the high-resolution data and at 3a" for the low-resolution data. 
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Fig. 5. — Position- velocity diagram taken along the major axis of the full-resolution data. The open squares 
show the rotation curve determined for the N and S side independently, while the crosses indicate the ve- 
locities for the SW extention of the galaxy. Contour levels are -0.135, 0.135, 0.27, 0.54, ... mJy beam'^ 




Fig. 6. — Vertical thickness of the H I emission expressed as the first moment of the z distribution at each 
position in the datacube with respect to the plane of the galaxy. Contour levels are 370 (white), 465, 560 
and 650 pc. 
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DGC 1807 




Right Ascension (J2000) Sight Ascension (J3000J 

Fig. 7.— Wide field total H I map with NGC 891 and the companion UGC 1807 (top left); Contours showing 
the total H I emission of the companion, contour levels 1, 2, and 5 x 10^" cm~^ (bottom left); velocity field 
of UGC 1807 (bottom right); p-Vplot of UGC 1807 taken along the kinematical major axis, contours 0.135, 
0.27 mJy beam"'. 
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Fig. 8. — Radio continuum contours overlaid on blue optical image taken from the Digitized Sky Survey. 
Contour levels are 40, 80, 160, 320, ... //Jy beam-i. 
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Continuum z distribution 




0.5 1 1.5 2 2.5 3 

z (arcmin) 



Fig. 9. — Normalised vertical density distributions of the radio continuum halo. The spatial resolution of 
the observations is indicated. The density distribution of the Ha emission as determined by Rand & Heald 
(priv. comm.) is indicated. 
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Fig. 10. — H I radial surface density for the thin disk of NGC 891 (squares) overlaid with the fitting formula 
in eq. 1. In the inner parts a compact exponential component has been added to model the H I ring. 
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Fig. 11. — Comparison between five representative channel maps for NGC891 (rightmost column) and 
models. From left to right the models are: 1) Thin disk; 2) Strong warp along the line of sight; 3) Disk 
flare; 4) Disk + co-rotating halo. Heliocentric velocities (km sec~^) are shown on top left comers of the data 
channel maps (vgys = 528 km sec~^). 
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Fig. 12. — Gas distribution in the region of the halo. Left panel: deprojected radial distribution at two 
distances (z = 2.8 and z - 5.6 kpc) from the plane (squares). The lines show the results of the fit with eq. 3. 
Right panel: scaleheight of the gaseous halo fthaio(^) (squares) with a power law fit (solid line). 
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Fig. 13. — Comparison between five representative channel maps for NGC891 (rightmost column) and 
models. All the models have a two component (disk+halo) structure with the halo lagging in rotation 
velocity with respect to the disk. From left to right they are: 1) Lagging halo with constant gradient; 2) 
Lagging halo with high velocity dispersion; 3) Lagging halo with a radial inflow motion; 4) Lagging halo 
with velocity gradient increasing in the inner parts. Heliocentric velocities (km sec"^) are shown in the top 
left comers of the data channel maps (vgys = 528 km sec~^). 
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Fig. 14. — Comparison between two position-velocity cuts perpendicular to the plane of NGC 89 1 (right- 
most column) and the models. The top row is taken 2.7' (7.5 kpc) North of the centre, the bottom row 1.0 ' 
(2.8 kpc) The models are the same as in Fig. [13] 
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Fig. 15. — Rotation curves as a function of z used in the modelling. In the inner region they become 
shallower for increasing z-distances from the plane (1.5, 3.0, 4.5 kpc). 
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Fig. 16. — Six representative channel maps for NGC891. The arrows indicate the filament (upper panels), two 
isolated "counter-rotating" clouds (middle panels), and the extension of the disk in the SW side (lower left panel). The 
lower right panel shows a channel map at high rotational velocity on the receding side. Note that both the halo gas and 
the extension in the disk are absent at these velocities. Contour levels are -0.4, -0.18, 0.18, 0.4, 1, 2, 4, 10, 20, 40 
mJy/beam. 



